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Abstract 



We study the two field model for reheating based on the scalar inflaton field, 
(p, and its interaction with another scalar field, through a Lagrangian term 
^g 2 X 2 93 2 ,much investigated for parametric resonance or preheating. Atten- 
tion is particularly on the quantum excitations of the inflaton field and the 
metric with a smooth transition from quantum to classical stochastic states, 
and with reheating followed right through from a specific inflation model to 
', a state including a relativistic (radiation) fluid. The excitations of the metric 

(but not the inflaton or x) are treated perturbatively and the validity of this 
approach is assessed. For this model our work points to the possibility of 
C (the curvature associated parameter for observed cosmic microwave back- 
q^. ground radiation anisotropies) changing significantly during reheating with 

QJ ', parametric resonance. 



I. INTRODUCTION 

Cosmic inflation theory combined with the theory of quantum fluctuations of the infla- 
ton field has provided a theoretical basis for the observed cosmic microwave background 
fluctuations (CMBRF) |JJ. In its simplest form with just one scalar field ip, the inflaton, the 
wave lengths of the relevant perturbations become much greater than the horizon distance 
during inflation and remain thereafter outside the horizon, uninfluenced by cosmic details, 
until they re-enter the horizon early in the matter era as density and matter perturbations. 
Among other things this neatly bypasses the complications of reheating when the inflaton 
field is assumed to thermalize into a radiation fluid. A parameter which is often used to 
calculate the CMBRF, for the relevant small wave numbers k, is the metric curvature com- 
ponent TZk which for these small k is approximately given by £ |2|,[| (with ( w —TZ); in 
the simple theory, having adiabatic dynamics of the physical quantities, this is constant 
throughout the rest of inflation, reheating and the radiation era. 

It has long been recognised that that the constancy of ( can alter in the presence of 
another scalar field, x- The question as to whether or not this occurs has to be addressed 
for the details of each particular theory and there has been considerable interest [B-Kj in the 
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effect of parametric resonance |9Hl5| on this question, and on the CMBRF generally. Quite 
a number of variations on this theme have been studied [12jl^,16jq| and some may well be 
significant for future directions. But basic scenarios are provided by the most investigated 
model where the interaction of the two fields is given by \g 2 X 2 ^ 2 'i parametric resonance 
may or may not occur depending on the value of the constant g 2 and other parameters of 
the theory. It is an object of this paper to investigate the above questions for the g 2 x 2( P 2 
theories where x has no initial classical component. 

A problem arising is treatment of the decoherence of the original quantum perturbations 
and the consequent stochastic classical variables |I7]]. In most treatments the CMBRF are 
evaluated with the matter era density and metric fluctuations (arising originally from the 
inflationary era) as stochastic variations [] 

In considering non-adiabaticity - whether or not accompanying parametric resonance - 
during reheating it is appropriate to consider wave numbers of up to the order of the Hubble 
parameter. All waves of x an d and tp% arise originally as quantum perturbations, where 
(f = (f + (fi with (f being the classical inflaton field. Then, treating all waves similarly 
they eventually smoothly decohere, in the way shown by Polarski and Starobinsky ||17|| , 
giving waves expressed by classical stochastic variables e(k) different for each value of the 
wave number k; the stage of decoherence depends on the number of the field quanta of a 
given k as a funtion of time. We shall retain the stochasticity throughout reheating and 
shall develop equations by taking averages over stochastic variables. Typically this will give 
rise to terms such as (x 2 ) where (...) denotes the average; this is the equivalent of taking 
vacuum expectation values in the quantum case, giving rise to a loop integral and being 
sometimes in this context called the Hartree approximation. 

We have used stochastic variables, with the ensemble averaging just mentioned; they have 
a smooth and well-defined descent from the quantum operators |17| . A simple alternative, 
which has been often used, is to insert at some time classical (non-stochastic) fields for x and 
(pi. Besides its intuitive nature that procedure has also the disadvantage of not being very 
well defined; so for the present paper we have adopted this alternative excursion into a more 
well defined treatment which has however its own approximation in the above averaging.^ 

In our treatment there is no need for a perturbative expansion in x an d <fi- In their 
occurence in the matter tensor Tj? of the Einstein equations, or in the scalar equations, the 
non-linear terms pose no difficulties since we have the averaging procedure and we get a set 
of simultaneous differential equations in the mode functions, as described in section III. For 
example {(pf) quickly becomes comparable with tp\ but this does not invalidate the equations. 
On the other hand we find it important to take account of the metric perturbation, ip, and 
this indeed we have to treat perturbatively and maintaining the perturbative validity imposes 



1 Such fluctuations have also been derived, by-passing decoherence, by keeping the Heisenberg 
representation quantum operators throughout time, and evaluating the final correlation functions 
in terms of the mode functions by taking (Heisenberg) vacuum expection values p^-pO[. 



2 The method of Khlebnikov and Tkachev [14| which is also based on stochastic fields is discussed 
below in section III. 
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a strong limitation on the parameters of the theory. 

To achieve a state where the x an d <p\ field densities compare with the faded classical 
inflaton field density is not to complete reheating. The usual view is that reheating is 
completed when the universe is dominated by a hydrodynamical fluid composed of relativistic 
particles. We go to such an era through the simplest mechanism of friction or dissipative 
terms with decay constants T x , r v whose contribution to the \i V 9 equations of motion 
respectively are proportional to T x x', T v ip' which operate effectively with rapid oscillations 
of the fields. Since oscillations, for example those associated with parametric resonance, are 
prominent features of the reheating process we expect such mechanisms or any dissipative 
decays, to have a significant influence throughout reheating. Such dissipative processes are 
entropic. 

These features provide a model smooth physical connection of the reheat period to a later 
era with a radiation fluid. To embed the reheat period precisely in a model cosmic history we 
need to specify the inflationary era to provide initial conditions for the reheating. We choose 
power-law inflation (partly for the convenience of an analytic solution) and this provides 
initial conditions both for the classical inflaton field tpo and its quantum perturbations <pi 
with the associated metric perturbation. During the inflationary period the \ field equation 
of motion has an important term g 2 X l Po- For any inflationary era in which (p is large, 
0(mp/ anc fc),and non-oscillatory this term ensures extreme suppression of any initial \ field. 
Thus the \ field, or equivalently the number of % particles, entering the reheat period is 
nearly zero. A smooth junction with the reheat equations of motion then completes the 
model embedding of the reheating in the assumed cosmic history. 

There now follows Section II in which we treat the relevant relations between quantum 
states and classical stochastic states. Sections III and IV give the reheat equations, their 
initial conditions arising from the inflationary model and some details of their solution. In 
Section V we discuss the results and we conclude with a general summary in Section VI. 



II. STOCHASTIC FUNCTIONS IN REHEATING 



Observational evidence [gj] on the cosmic microwave background fluctuations tends to 
support the theory that these are due to stochastic perturbations of the cosmic density and 
metric arising from quantum perturbations of the early classical inflaton field, fo(r) [|TJ. 
These can be written as 



^i( x > r )= / — -r r[c(k)^ k (r)exp(2k.x) + h.c] (1) 

J (27T I 2 



</?i(x, r) is an operator in the Heisenberg picture, the time dependence being wholly con- 
tained in the mode functions ^{t), and c(k) is a time-independent quantum annihilation 
operator such that [c(k), c(k')t] = 5 3 (k — k'). ^fc(r) is a mode such that its r dependence, 
for r — > oo, is given by 

(fk(T) oc a -1 exp — ikr (2) 

We have adopted a conformal time metric whose unperturbed form, having scale factor 
a = a(r), is 
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ds 2 = —a 2 dr 2 + a 2 5ijdx l dx J . (3) 

The most studied form of interaction has been with another scalar field x with the 
interaction potential given by 

v in t{v,x) = lgVx 2 (4) 

where (p = ip + {pi and x is i n origin a quantum field. We adopt this as our representative 
of the basic scenario. Like most other authors in this context we do not consider the case 
where x has an aboriginal classical scalar field part; this would change the model to be more 
sophisticated like some variety of two-field inflation. Thus we write % similarly to (pi as 

/d^ k 
—— s-[6(k)Xk(r) exp(zk.x) + h.c] (5) 
(27TJ 2 

where 6(k) is a quantum annihilation operator , [6(k), &(k')t] = <5 3 (k — k'), and Xk( T ) is a 
mode function such that its r, k dependence, for r — > oo, is given by 

Xk( r ) oc a -1 exp — ikr/\^2k (6) 

Thus we regard x as originating like ipi. In the absence of a deeper theory, and in view of 
the seeming success in CMBR calculations of the theory of (p\, this seems a most reasonable 
ansatz. The <^fc(r) develop through the inflationary era and so do the Xfc( r ) mode functions. 

Now it is quite legitimate, for linear coupled equations of motion in the various quantum 
fields, to carry this development right through the reheat and radiation eras to the origin of 
the CMBR while evaluating the Heisenberg operators and finally taking vacuum expectation 



values of correlation functions |2^ , p!9|j20|| . In dealing with non-linear coupled equations it 
can be appropriate to make the transition to the more usual classical stochastic picture of 
the fields. Let us discuss this using as an example the scalar field equation for x which is 

X (x, t)" + 2(a'/a)x(x, r)' - V 2 X (x, r) + a 2 [M 2 + g 2 ^) + ¥>i(x, r)) 2 ] X (x, r) = (7) 

This is, from Eqs.(|l]||), a non-linear equation in the quantum operators where the non- 
linearity arises from the g 2 coupling of the operators in (p% with those in x- We eliminate 
the great quantum field theory difficulties thus arising by taking the Hartree approximation 
- discussed further below - on the c, in ip\. That is we replace (<^o(' r ) + ¥i( x i T )) 2 by its 
vacuum expection value (((Po(t) + ipi(x,r)) 2 ) = (po(r) 2 + (<^i(x,r) 2 ) 

(i Pl (x,r) 2 } = (2n)- 3 [d 3 ^^ (8) 



In field theoretic terms this is a loop integral and as stated in Eq.(|18D below the same 
integral form is maintained when we go from quantum operators to stochastic variables. Now 
take the Fourier transform of the resulting approximation to Eq.(|7]), leading to an equation 
linear in the creation and annihilation quantum operators bk and b k . The coefficients of bk 
and b\ must each vanish giving two equations which are complex conjugates. We could have 
reversed the order of the derivation and taken the vev after the selection of the coefficients. 
The resulting x mode equation is 
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Xl + 2{a'/a)x' k + (k 2 + a 2 M 2 ) Xk = 



(9) 



M 2 (r) = M 2 + <rWr) 2 + ^(x, r)) 2 ) (10) 

It is important to note that \ is complex, as is the corresponding (pi in its equation, and the 
imaginary and the real parts are not constant multiples of each other. We have spelt out 
this derivation in order to point up the similarities and the contrast to the semi-classical, 
stochastic, case which now follows. 

Polarski and Starobinsky |17 have studied the transition from the quantum to the semi- 



classical case. For the Heisenberg picture which we have they find that when the mode 
functions grow large enough so that Planck's constant, h, can be neglected then the creation 
and annihilation operators may be replaced with appropriate Gaussian variables to form an 
equivalent stochastic field. In our formalism Eqs. (|I],|5]) are replaced by 

/d 3 k 
—— T [e(k)^ k (r) exp(zk.x) + e*(k)y£(r) exp(-ik.x)] (11) 
(2n)2 

X (x,r) = / — - T [d(k) Xk (r)exp(ikx) + d*(k) X *(r)exp(-ik.x)] (12) 

e(k) and d(k) are time-independent separately 5-correlated Gaussian variables such that, 
where (....) denotes the average, 

(e(k) e *(k')> = (d(k)«T(kO> = ^ 3 (k - k'); (13) 

(e(k)e(k')) = (d(k)d(kO) = 0; (14) 

(e(k)rf(k')) = (e(k)d*(k')) = 0, (15) 

that is the e variables and the d variables are uncorrelated. 

If now we insert Eqs. (|11|,|1^) into Eq. ([?]), take the Fourier transform, and average over 
the Gaussian variables e we obtain, analogously to Eq.(U), 

^ + 2(a'/a)x k + (k 2 + a 2 M 2 )x k = (16) 

Xk = rf(k) Xk (r) + d*(k) Xk (r) (17) 
where M 2 is given by Eq.([T0|) since 

(Cfl) = (27T)- 3 f d 3 kif k ifl = (cpDo (18) 



The coefficients of the uncorrelated variables (d—d*) and (d+d*) must each vanish separately 
and we can thus obtain 
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Xi: + 2(a'/a)xL + (* 2 + a 2 M 2 ) X ^ = 



(19) 



where Xk is complex and because of Eq.QTSD this is the same as Eq.(|]).Thus we have a 
smooth transition from quantum to semi-classical or stochastic. 

Polarski and Starobinsky have also shown that, in the semi-clasical regime, the mode 
functions can be made real by a time independent phase transformation and they formulate 
the theory with somewhat different Gaussian variables e from those of Eqs.( |l4l , |l5|) . We can 
transform our classical stochastic formulation to that of Polarski and Starobinsky: 

We can write Eq . (fTT|) as 

d 3 k 
(2vr) 



— - r[e(k)^ k (r) + e*(-k)^_ k (r)] exp(zk.x) (20) 

I y.TT I 2 



and we can time-independently transform our variables so that <pk becomes real .Then using 
the fact that (p^ = <£>_k we get 

— - T e(k)(^ k (r)exp(ik.x) (21) 

(27TJ 2 

where e(k) = e(k) + e*(— k). Thus 

e(k) = e*(-k) (22) 



(e(k)e*(k')> = 5 3 (k-k') (23) 

where the last equation is derived using Eq.(|T3|). Eqs. fl2"I| - |23D are those of Polarski and 
Starobinsky (except that we have the factor a~ l (r) incorporated in our mode functions). 

To summarise: using the Hartree approximation the scalar quantum and classical 
stochastic mode equations are the same, with complex mode functions, so there is a seam- 
less join. In the classical stochastic regime a conversion to real mode functions can be made 
yielding completely real equations, but these equations are invalid in the quantum regime. 
For this reason we prefer to use complex mode functions in reheat studies as not all modes 
may be classical stochastic through all the relevant period. 

We note that if £{k),T)(k) are semi-classical mode functions, which must be of dimen- 
sion (mass)~2, then we cannot write equations such as + .... + g 2 J £(k — k')r)(k' — 
k")r](k")d 3 k'd 3 k" = which if valid would contain different and interesting rescattering ef- 
fects. Apart from being dimensionally inconsistent we would have had to take averages 
over the stochastic variables in their derivation and this, as previously above, would have 
eliminated all such terms. Of course such equations are perfectly possible if £, rj are fully 
classical mode functions of dimension {mass)~ 2 and indeed have often been written down in 
the resonant reheating context. But to obtain such such equations, starting from a quantum 
origin of the modes as we do, requires specific transformations and we are not aware of 
any suggestions for these for k ^ 0. Naturally if = d(k)xkiV(k) = c (k)<Pk then such 
equations are valid, but the physical solution involves solving a statistical ensemble of equa- 
tions, got by sampling d, e and finally averaging - a seemingly very lengthy procedure. The 
method of Khlebnikov and Tkachev |TJ|], mentioned also below, seems to provide at least a 
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partial and more achievable way of implementing the stochasticity in the semi-classical case 
without resort to the Hartree approximation. However having a system of 7 simultaneous 
equations with complex mode functions we shall implement the quantum or semi-classical 
equations which we have illustrated above. 

As a final comment we note that it is dimensionally obviously inconsistent to write such 
coupled classical equations as the above and loop integrals such as ([18]) with the same mode 
functions. 



III. THE REHEATING EQUATIONS 

We now formulate the equations of motion for the reheat period. The Lagrangian for 
the scalar fields is 

tfxyf^p^ + \ X a X, a ~ V{y) - V{ X ) - V int {cp, x)) (24) 

where <p = <po + <pi and Vi nt is given by Eq.(f4j). 

In the perturbed metric, g a ^, we use the longitudinal gauge for the perturbation 

ds 2 = a(r) 2 (l + 2(j))dT 2 - a(r) 2 (l - 2^)5 ij dx i dx j (25) 

Unusually, the scale factor is calculated using the full reheating dynamics as below, Eq. ([33]) . 
Unlike other treatments, excepting Bassett et al. ||, we include the metric perturbations 
as dynamical variables in the coupled equations of motion. In this physical system, includ- 
ing the hydrodynamical variables introduced below, there are no space-space non-diagonal 
components in the energy-momentum tensor. Hence <p — ip P3 and we denote this metric 
variable by ip. 

We also wish to take acount of actual reheating in which the <p and the produced x 



particles, which may be massive [^4], are replaced by a thermal gas of relativistic particles, 
giving the radiation era. We make no attempt to imagine details of this process but adopt 
the simple friction mechanism which puts all these details into a black box. The use of 
this mechanism in general cicumstances has been criticized, it being claimed that it is only 



appropriate when the decaying field is undergoing fast oscillations [13]. In the calculations 
that we shall describe it is in these circumstances that the main production of a relativistic 
gas occurs. We describe this gas by the usual hydrodynamical variables p, the density, and 
p the pressure and in the equations of motion which follow we take that equation of state, 
p = p/3, which is appropriate for relativistic particles. 

We divide the hydrodynamical density into two components 

p(x, r) = po (r) + pi (x, r) (26) 

where poij) is the homogeneous background density and pi(x, r) is the inhomogeneous part 
of it; the notation 5p = pi/po is the usual usage. As we shall see in the equations of motion 
pi mainly arises from the inhomogeneous part ipi of the inflaton field in accord with the 
usual theory of the cosmic microwave background fluctuations. 
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We shall now formulate the equations of motion, being those Einstein equations and 
scalar wave equations arising from the Lagrangian, Eq.(24), and metric, Eq . (p5|) , with the 
addition of the hydrodynamical terms. The energy-momentum tensor is 



T? = y^yv + x ^x,u - [-^.a + -X' a X,« + V(<p) + V( X ) + V mt - p]<5£ + (p + p)u»u u 

(27) 

where is the fluid 4-velocity. 

First we write 3 spatially homogeneous ,or background, equations taking the average 
over the Gaussian variables as in Eqs (|l3l - |l8l) . We use the notation ' for djdr and we take 
the individual field potentials in reheat to be: 

V^) = \m\ 2 -V[x) = \M\ 2 . (28) 

After Eq.([37]) below it is noted that tp% and if) have the same stochastic variables. Thus 
defining the density and pressure homogeneous parts of the energy-momentum tensor as 
Pt(t) = —(Tp) and Pt{j)8] = (Tj) we find, after ensemble averaging, that 

Pt(t) = ^[V + a 2 m 2 (if 2 + (<p 2 )) + a 2 M 2 ( X 2 ) + (<fi&) + (xj) - Vo<^>] + Po, (29) 



Pt(t) = ^ 2 iv- a 2 m 2 {v 2 + <tf» - a 2 M 2 ( X 2 ) - (^.»/3 - (x? 4 )/3 - VofoWl + A>A 

(30) 

»7 = + K) + (X' 2 ) (31) 

rh 2 = m 2 + g 2 (x 2 ) . (32) 
and we have the Friedmann equation 

(a /a) 2 = ^a 2 PT (r). (33) 

As in Eq.fll8|) the averages are independent of x ensuring the spatially homogeneity of 
the R.H.S. of Eq.(^). We also see from Eq.([Tj|) that Eq.([33|) is equally valid in the quantum 
regime since the formalism ensures that the vev of the quantum operators equals the average 
over the stochastic variables. 

The other two spatially homogeneous equations are for (po and po- Here we introduce the 
arbitrary decay constants and T x which serve to create the hydrodynamical radiation 
gas specified by p — po + pi,p — p/3. In Appendix A we show how these are introduced in 
a way consistent with the Bianchi identities and how the equations for p , pi, <p , <pi, and x 
are subsequently derived. 
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<p'o + 2(a'/a)<p' + a 2 rh 2 p - 4<VVi> - 4(^VVi> + 2a 2 m 2 (^i> = -a^M + 2(M)) (34) 

We note the combined reaction of the inhomogeneous fields tpx and ip on the homogeneous 
inflaton field through the ensemble averaging, similarly to that in Eqs. ( p9|) and (|30|). 

p(, + 4(o7o)A) = a- x r v (^ 2 + K» + a-^^x' 2 } - 2a~ 2 ^V 2 ^) + 4(^'pi> (35) 

There remain the 4 spatially non-homogeneous equations which we write in the /c-component 
form. For \ and (pi these components are specified as the complex mode functions Xk and 
p k of Eqs. (EL [IT] , ) . Their wave number dependence, given by the succeeding equations, is 
only on k = |k| . 

Xl + 2(a'/a) X ' k + (k 2 + a 2 M 2 ) Xk = -aV xX ' k (36) 
where M 2 is a function of r given by Eqs. (|l0| , [T8| ). 

(fl + 2(a'/a)(p' k + (k 2 + a 2 m 2 ) Vk - 4(p' ip' k + 2a 2 m 2 V ^ k = -aT v ((p' k + 2(p' Q ip k ) (37) 

Though not indicated these equations actually hold for each separate value of k; thus con- 
sistency of Eq.(^) indicates that ip^, the mode function of the metric perturbation should 
be associated with the the same Gaussian operators (or, in the quantum regime, with the 
same quantum operators) as <pk in Eq.flTTp. Thus 

</>(x,r) = I -^ 5 -[e(k)^ k (r)exp(ik.x) + e*(k)^(r) exp(-zk.x)] (38) 
J (2tt)2 

where since t/>(x, r) is dimensionless ip^ has dimension (mass)~ a . This associates the metric 
perturbation with the inhomogeneous part of the inflaton field without assigning priority to 
either. But the stochastic variables of Xk are independent. Thus no terms in ifa appear in 
Eq. fl36"D ; they are forbidden through ensemble averaging.. 
The mode equations for ip and p\ are 

M + (3^ + (a'/a)ip k )(a'/a) = ^a 2 (p T + 3p T )^ k + ^Ga 2 b Vk (39) 

where 

8Pk = \[-£if>k + <Po<Pk - 2<p' k (^<p' 1 ) - a 2 m 2 (p <p k } + p k /3 (40) 



with 



e = p' 2 + K 2 > + (x' 2 ) + (<Pim,i)/3 + (x,<X,i>/3. (41) 



p' k + A{a'/a)p k - 4poVl = -^[(^G)- 1 ^ + (a'/a)^ k ) - (p' + 2<Vi>M + 21>^' fc 



(42) 
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Eqs.(|39D-(f42|) can be deduced by extending the space-time and space-space gauge invariant 
Einstein equations given by Mukhanov et al. |23|1 ; the term in square brackets in Eq.([42|) 
arises from the validity of the space-time equation and vanishes if po — 0. 

From the structure of the above equations the Fourier transform of p\ (x, t) has the same 
Gaussian variables, e(k), as have <pi(x, t) and ^(x, t), having dimension (mass)i] the 
quantum operators and consequent Gaussian variables of the \ field are independent of the 
others. This point also has consequences when we use vacuum expectation values or averages 
over the stochastic variables as noted above in connection with Eq . (j36[) . 

A number of points arise on our equations of motion in the context of the large body 
of work on parametric resonance in reheating since the subject was introduced by Traschen 
and Brandenberger M: 

i. Scale factor. We find the scale factor a(r) by Eq.(^) this being the appropriate 
calculation which includes all the relevant dynamical variables, and thus the influence of all 
the scalar fields. Generally an approximate ansatz for a has been used in other paper, an 



exception being the work of Boyanovsky et al. JT^|. Our equation makes a difference to the 
form of a but this change does not seem to have a critical influence. 

ii. The loop integrals (Hartree approximation). For quadratic field forms, in the classical 
stochastic regime, we have used the Gaussian average which is a smooth continuation from 
the quantum field theory vacuum expectation value. This has been called the 'Hartree 
approximation' and has been much used in parametric resonance work. (We give some details 
of our evaluation in Appendix B.) We note that, for example, the term g 2 fl(x 2 ) corresponds 
to a first order single loop calculation in quantum field theory as discussed by Kofman et al. 
T2|| . They have compared the magnitude of this with that of an approximate evaluation of 



the single loop second order in g 2 term, and have concluded that this could be of comparable 
magnitude for some resonance modes they studied. In common with Kofman et al. and 
other authors we have not included this difficult, though maybe significant, correction. To 
include suchlike corrections we would have to add the results of the appropriate field theory 
calculations to Eq. (p4|) , thus manufacturing an effective Lagrangian. 

As discussed in the last paragraph of section [TJ our formalism does not include rescat- 
tering corrections of the form / C,(k — k')i](k' — k")i](k")d 3 k'd 3 k" as our regime is stochastic. 
Other authors have written down, though perhaps not fully implemented , equations 

with such terms which seem to require the fields to be already fully classical, in the sense 
noted in section || together with a question on their provenance. We have found that our 
emphasis on the Hartree approximation is closely related to the viewpoint of Boyanovsky et 
al. [TJ. 

Hi. The metric perturbation. The coupled equations include the metric perturbation, ip. 
As noted above its quantum or stochastic variables are the same as those of (p± whose equa- 
tion of motion (^) has an important coupling to ip arising from the metric interaction with 
y?o- The equation of motion fl39|) of ip with Eqs.( ^0| , p] ) show its coupling to all the other 
fields. In the derivation of the coupled equations ip is treated as a small perturbation to 
the metric, this being in contrast to (f and x which are not treated perturbatively. So the 
validity of this has to be maintained throughout reheating and we find this gives significant 
restrictions on the parameters of the theory. Much work involving solving equations of mo- 
tion in this context has ignored the metric perturbation but Bassett et al. H have discussed 
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and emphasized its importance. 

iv. Stochasticity. Khlebnikov and Tkachev [14], as mentioned also in section |J have 



an entirely different approach to the equations of motion, which they solve on an x-space 
lattice. Their method takes account of the stochastic nature of the variables, as exhibited 
for example in Eg. (|TTD , by taking appropriate initial conditions. For each discretized k- 
value they take a random value (chosen from an appropriate distribution) of for example 
the variable e(k)<£>k(r) °f Eq.fllT]). The initial value of y?(x, r) is then got by using the 
discretized Eq. (|TT|) to sum over all those random values. Thus though only one value is used 
for each k each x contains a large random sample. After solution in x-space the variables 
can be re-Fourierized to get the k components. Because of the non-linearity of the x-space 
equations this mixes the initial k components and is not equivalent to solving the k equations 
directly. We note that the fields must already have passed through the quantum stage into 
the classical stochastic phase; the authors do not include the metric perturbation . 



IV. SOLVING THE EQUATIONS 

A. The Microwave Background Fluctuations 

The present inflationary theory of the origins of the cosmic microwave background fluc- 
tuations has so far been successful. In its simplest form a perturbation, S(fik(t) of quantum 
origin, in the inflaton field developes in the inflationary era in association with the curvature 
perturbation, TZk(t). At a certain time, t*, in the inflationary era before reheat begins, the 
modes of relevance are such that k/a becomes equal to, and then rapidly less than, a/a = H. 
Since then \hysicai = 2ira/k ^> H^ 1 the latter being, in this context, the Hubble radius and 
this is sometimes expressed as the perturbations having passed outside the horizon. At 
around the beginning of the period of matter domination if -1 again becomes greater than 
^physical- T^k(t*) = — [H 5 (p k / (p]t=u and, with certain assumptions in the simple model such 
as adiabaticity, IZk remains the same during that whole period when X p h ys icai 3> H~ x and 
equivalent || to the parameter ( 0.^3|.||.(2i| : 



( k = ^(H-^ k + $ fc )/(l + w) + $ fc = -TZ k (43) 

where <Dfc is the gauge invariant metric perturbation (equal to the longitudinal gauge metric 
perturbation ip k of Eq.(|25|)) and w = p/p being the total pressure to density ratio. This 
scenario relates the CMBRF rather directly to the inflationary period bypassing reheating 
complications. 

Now the model which we are studying has an extra scalar field, with an interaction g 2 x 2( P 2 
known to be capable of giving rise to parametric resonance (preheating [|n],|l2|]) in reheating 
. We consider how far and under what conditions, the constancy of ( for \ p h ys icai H^ 1 
remains valid. For this purpose we investigate the conclusion of reheating by particle decay, 
as mediated by the and parameters and for comparison we also look at cases where 
these are zero. 

In what follows we shall largely be engaged with modes which have wavelengths which do 
not become larger than if -1 before reheat. So we should have an estimate of the magnitude 
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of all relevant modes at the beginning of reheat and this will be dealt with in the following 
subsection. 

B. Initial Conditions at Reheat 

1. x- field 

Whatever we may postulate in detail for the single field potentials V((p),V(x) through 
the inflationary and reheat periods, it is a hypothesis of our model that the g 2 f 2 x 2 is the 
significant interaction term in the Lagrangian. This has some immediate consequences for 
the magnitude of x at the (fuzzy) interface between inflation and reheat. A simplified version 
of the mode equation, Eq.fl3"ED with Eq.([H]), in the inflationary era is 

(a X k)" + [k 2 + a 2 M 2 + g 2 a\l - a"/a](a X k) = (44) 

and ipl is greater than, or of the order of, m 2 Pl through most of the inflationary era. Thus 
the term in square brackets is large and positive resulting in a quasi-periodic type solution 
for a^fc, which has a primordial value given by Eq.([|). The many efold increase of a during 
the inflationary era indicates an exceedingly small value for Xk a t the beginning of reheat 
p6| , |27|| , which may be as little as 10~ 50 of its initial value. This is an important qualitative 
feature. Firstly for such small mode functions in the beginning of reheat the transition 
to classical stochastic functions cannot yet be made and the quantum complex formalism 
should be retained. We have taken the initial value of x to be W~ n m pl ^ nck where n is of the 
order of 30. 

2. Background metric and the classical inflaton field 

For the reasons just outlined we shall assume that it is valid to neglect any influence 
of the x-field during inflation on the classical background field ipo and on the quantum or 
stochastic part ip\. So we shall work with the standard single field inflation formalism to 
formulate conditions at the beginning of reheat for both ipo and ip\. It will prove convenient, 
as we shall see in 4 below, to have a potential, V(ip), of different form in the inflation 
era from that in the reheat era. In a perfect physical model these should be joined by a 
smooth transition form but as we show we can substitute a sudden (phase change) transition 
for the purpose of finding the initial conditions of reheat. First we deal with the classical 
background field which we denote just by ip in the rest of this section. 

At the boundary between inflation and reheat we can start the reheat era with the same 
values of <p and ip' as those at the end of inflation and we can likewise impose continuity of 
the potential V (for an example see 4 below); of course we also take a continuous. It follows 
that a' is continuous since 

a 2 = {a>/a) 2 = ^[\ v ' 2 + a 2 V} (45) 
We can show that a' is continuous as follows. 
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2aa' = — fo>V" + a 2 ip'(dV/d<p) + 2aa 2 V] 



(46) 



Since 

tp" + 2atp' + a 2 dV/dip = (47) 

it follows that 

7 = 1-4 = 47rGV/a) 2 (48) 
or 

where 7, thus defined, is a quantity we shall need in the next section. 

Thus from the imposed conditions at the beginning of reheat it follows that a', or equiv- 
alently a", is continuous at the boundary. Thus so is the background curvature, in addition 
to the metric. 



3. (pi and the metric perturbation 



Having a satisfactory background continuation we can now deal with the continuation of 



the quantum perturbations of the inflationary era. Deruelle and Mukhanov |28| have used 
the Lichnerowitz curvature conditions to give the necessary continuity conditions on the 

ip, the metric perturbation 



metric perturbation. These can be summarised as follows f20|j2 
and T should be continuous across the boundary where 



T = (ip'/a + V - V 2 V/3a 2 )/7. 



(49) 



Since 7 has just been found to be continuous the last condition can be reduced to the 



continuity of ip' + aip, which in turn using the time-space Einstein equation [23 



ip' + aip = 4:7iG(p' ip 1 



(50) 



can be reduced to the continuity of cpx. 

Thus from the Lichnerowitz conditions one finds that the values of the mode functions 
for the metric perturbation ip and for the non-homogeneous part, (pi, of the inflaton field at 
the beginning of reheating after the phase change are the same as those at the end of the 
inflation era; ip', given by Eq.([50|), is also continuous. <p[ does not need to be continuous 
and in fact is not. Its value can be found in terms of the continuous fields by using the 
time-time Einstein equation at the beginning of reheating (which differs from that at the 
end of inflation). At this very beginning of reheat, as at the end of inflation, we validly 
use linear perturbation theory in ipi and we can check that the perturbations are adiabatic. 
This means that we expect ( to be constant in the first part of the reheat period, which is 
what we find as discussed below in section IV C. 



The values of the above quantities with cp , cp' and a = (a' /a), plus the input of the very 
small X-, x' 1 provide all that is needed for the initiation of the reheat equations. 
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4- Input from power-law inflation 



Rather than making some plausible input data from a generality of of inflationary cal- 
culations we have chosen to use a specific inflationary model. This is power-law inflation 
with an exponential potential: V = Uexp(-Xip) where U is a constant. In this there are 2 
parameters free to be chosen. These are the specific power of the inflationary scale factor 
a, and the magnitude of the classical scalar field <po at the junction of inflation and reheat, 
which we define as where we begin to use the reheat equations of section Til. 

An advantage of this power-law inflation is that the solutions are analytically expressible 
2T| . For the scale factor a and the inflaton field ipo 



a oc (Tj — t) p ,lp' = —a\/n, (51) 

where Tj > r(end inflation) is a constant and p = 2k/(X 2 — 2k). With X 2 / k <C 2, p is near 
to —1 and there is power-law inflation. The imposed continuity of V at the boundary gives 
Uexp(— \(po) = \m 2 ^\\ thus, given m, A, we can arbitrarily specify ipo at the boundary and 
this condition then merely specifies U. This triviality in making an arbitrary choice of ipo 
at the beginning of reheat is useful. 

Also, in the inflation era, the perturbation <^i(x,t) and the perturbation -0(x, r) to the 
metric can be expressed in terms of a quantum field function 

Z3^| Kk)/i k (r) exp(ik.x) + h.c] (52) 



whose k- modes satisfy 

fi'l + (k 2 - a"/a)/j k = (53) 



This function emerges in the synchronous gauge, a derivation due to Grischuk [|19|J20| ;one can 
convert the resulting formulae to our longitudinal gauge quantities by well known techniques 
2|,[2(|. These yield 



*Pk = a 1 V-^ a y/h\(tfk ~ a ^k)/k 2 (54) 

ip k = a-\2ky 1/2 [fi + a^fjf - afx)/k 2 ] (55) 

where 7 = 1 — a' /a 2 = (p + I) /p. Since a" /a = —^{^i — t) -2 , Eq.flSSD can be expressed as 
a Bessels equation in terms of 

k 

y = fc(r< - t) = — \p\ (56) 
an 

and the solutions which satisfy the asymptotic conditions for a scalar field as in Eq.(H) are 

jjiy 1 1 

Mv) = \j^{Jn- iY n ) exp[-i(-wr + ~tt)] (57) 
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where n — 1/2 — p. 

For k/a <C H (many times over fulfilled by the k relevant to the CMBRF observations) 
/ifc(y) oc y p oc k p and if this /c-dependence passes unaltered through the reheat era into 
the radiation era it yields a power-law spectrum which scales as k n ~ x where n, the spectral 
index, is given by n = 2p + 3. 

Besides (po and the power p there is one more parameter that determines the input from 
inflation and that is the m of the inflaton field potential in reheat, V((p) = \m 2 ip 2 . This is 
because it is necessary to know the value of y, defined by Eq.fl56|), in order to determine the 
size of the perturbations. From the conditions of continuity of section [IV B 2| it follows that 
at the junction of inflation and reheat 



H = V^rGrwp y/p/(2p - 1) (58) 
and one can use this in Eq.([56|) to give 

»=s?w (59) 

To summarize: rmpo and p determine the input of ip and (pi to the reheat era and the 
reheating equations of Section III determine subsequent developments. 



V. RESULTS 

Our interest is in the qualitative features of reheat in a model compatible with known 
observations rather than in making precise comparisons with data. So we are interested in 
values of the parameters that roughly supply the usual requirements of inflation and the 
magnitude of the CMBR fluctuations. 

Unless otherwise stated the results we quote are for the values m = 10 _7 ,y?o = 0.3, p = 
— 1.1, M/m = 0.02. Throughout we quote dimensionful results and parameters in units such 
that U — c — G — 1. 

The critical parameters are the coupling constant g 2 and the frictional decay constants 
T x and r^, of the x an d p>- The quartic coupling constant g 2 governs the number of x produced 
and in particular the parametric resonance. We have considered values in the range given 
by 20,000 > g/m > 1000. For smaller values there is no appreciable particle production. 
The range of consideration for the frictional decay constants was 0.1 > r x /m > and 
0.005 > T^/rn > 0. Larger values of T x tend to suppress parametric resonance while larger 
values of may be unrealistic in the context of the development of the classical inflaton 
field. We have by no means done a complete scan over this parameter range but in each 
particular case we have looked at (i) the existence or not of parametric resonance (ii) the 
development of the metric inhomogeneity ip and also ( (iii) the post reheating state; and 
other significant features. To illustrate significant features the figures show a case with 
parametric resonance in the production of \ particles and r x /m = 0.001, T^/m = 0.0004 
so that no x particles are finally left. If both T's are zero then reheating (as opposed to 
preheating) does not happen. But it is interesting to inspect this case (also shown in the 
figures) to look at what the possible variation of ( might be if the decay of the x an d f 
particles is significantly slower. 
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A. The Metric Perturbation,?/; 



As noted in section |ITT] the metric is given by ds 2 = a(r) 2 (l + 2tp(x, r))dT 2 — a(r) 2 (l — 
2tp(x, r))5ijdx l dx : ' .From the end of reheat tp develops through the radiation and matter eras 
so that the magnitude of the components tp(k, r) for very small k is a determinant of the 
microwave background perturbation. 



1. The magnitude of tp 

The validity of the perturbation formalism depends on tp(x, r) being small and for each 
case we should investigate this, noting that it is not a quite straightforward concept. 

We shall discuss tp as having developed into a classical stochastic field, Eq. (581) . The 



indeterminancy represented by the variables e forces us to consider the average over the 
product of these variables so that we evaluate < (tp(x, r) 2 > the result being 

(tP(x,r) 2 ) = (2vr)- 3 / d 3 kMr)Mr)*, (60) 



the same for every value of x. We require that y< (tp(x,r) 2 > be small compared with 
unity, since the relevant metric coefficient, Eq.^5), is a(r) 2 (l + 2tp). Our viewpoint is that in 
any particular case satisfaction of this requirement forms sufficient justification for the per- 
turbative approach, because the only way we can mount a comparison of the revised metric 
coefficient with a(r) 2 is when we consider the basic equations to be those in configuration 
space, and then indeed the revision is just a perturbation through all space-time. It is true 
that if, having written the equations of motion in configuration space, we then project into 
fc-space to solve them, as we do, then certain of the tp{k) will be very large. But at this 
stage we are operating translation into /c-space as a device to aid solution. Also then the 
comparators from a(r) will be much larger, these translating either into a fc-space ^-function 
or, in the non-ideal case, into a function very strongly peaked round the origin in fc-space. 

In the cases illustrated we see from Figure 3 that the requirement is very well satisfied 
throughout and from Figure 4 that the requirement is also satisfied though significantly 
better in the first part of reheating. 



B. Production of \ an d Inflaton Particles 

It is a generic feature of this model, and probably of most other similar models, that the 
inflaton particles or inhomogeneous field, (pi, that enter the reheating period from inflation 
attain an energy density comparable or greater than that of the classical inflaton field, 
(po, during the preheating period and thereafter. In their equation of motion, (57|), we note 
the bilinear couplings to ipo and the metric perturbation, tp. Thus in preheating ip± cannot 
be treated as a perturbation and we cannot neglect factors like {(pf). 

The occurence of parametric resonance for \ is as expected sensitive to the value of 
g, and it is also, naturally, sensitive to the value of T x which tends to inhibit it. In the 
case which we illustrate we see from Figure 1 that parametric resonance occurs markedly 
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but rather briefly and that afterwards the \ particles almost vanish under the influence of 
the friction decay into a radiation fluid. As previously noted x an d ip have independent 
stochastic variables. 



C. Magnitude of ( During Reheating 

It has been found useful to define the parameter ( of Eq.[J3] where Cfc ~ — T^hi the 
curvature perturbation 0^3|. f\ For adiabatic perturbations with k 2 <C a 2 H 2 , is constant. 
For example, in eras where there is only one scalar field or only one hydrodynamic fluid, then 
the perturbations are adiabatic. This holds before the reheating period because the scalar 
field x is negligable, as we have seen above. However in the reheating period we have two 
scalar fields (pi, x and also the radiation fluid (which arises from dissipative processes) giving 
rise to an entropy perturbation, [Q,and thus the theorem does not apply. Importantly, in 
addition, (pi can no longer be treated as a perturbation and we do not do so as has been 



emphasized in section [VB . 



In the equations of motion the wave number appears as k 2 and we evaluate ( for a value of 
k such that k 2 < lO~ A a 2 H 2 throughout the reheating period. We see from Figs. 3 and 5 that 
C has in the radiation era modestly decreased from its value at the end of the inflationary 
era. When we switch off the dissipative decay into a fluid, so that there is preheating without 
reheating, ( now increases markedly as shown in Figs. 4 and 6. We also see from Figs. 5 
and 6 that ( is initially constant as it should be since (pi is then small enough to be treated 
to first order only. 



VI. SUMMARY AND CONCLUSIONS 

We have studied a single field (<£>) model of inflation, with another scalar field (x) which is 
naturally quiescent during inflation and becomes active post inflation with the possibility of 
parametric resonance and preheating as the classical part, <po, of the inflaton field decreases 
and oscillates. We have included reheating through friction-type decay mechanisms of the 
scalar fields into a gas of unspecified relativistic particles. We have adopted a method having 
stochastic variables naturally succeeding the quantum operators of the early inflation era 
and have developed equations for the post-inflation period by taking averages over the 
stochastic bilinear forms in the scalar fields. This enables a non-perturbative treatment of 
the stochastic scalar fields x an d p>\ (where ip = ip + pi) though the metric perturbation 
field, t/j, has to be treated perturbatively. 



C being a linear function of the curvature perturbation if) arises with the same quantum operators 
as the non-homogeneous inflaton field pi as was noted for ip in section III. In the successor classical 
stochastic variables the mode functions of p>\ , ^ and £ are complex but as we have noted previously 
can be converted to real by a time independent phase transformation |jTfl1 . It is this real quantity 
which we quote here for 
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We find that the non-homogeneous scalar field (which might also be described as bearing 
inflaton particles), ip±, remains strong with an energy density which competes (according to 
the values of the decay parameters) with that of the gas of relativistic particles, during a 
large part of the reheating period. 

There are cases of the parameters in which the metric perturbation tp markedly resonates, 
at much the same time as a parametric resonance in x, and exceeds the perturbative limit. 
We have to reject such cases from consideration though we cannot say whether the failure 
be mathematical or physical. There are, however, many valid cases. 

The curvature parameter, (, for wavelengths relevant to the cosmic microwave back- 
ground fluctuations occurs very conveniently as a constant in non-entropic single field in- 
flation models. In this present more complicated model it changes from its value in the 
inflationary period by up to a factor of the order of 10 in the examples presented in this 
paper. 

We thank Andrew Liddle, David Lyth and Luis Mendes for valuable conversations. 
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FIGURES 




FIG. 1. Logarithmic plot of energy densities of x (A)and ip$ (B) versus the number of e-folds 
of expansion in the reheat era. The parameters are g/m = 10 4 ,r x /m = 0.01, Y^/m = 0.004, with 
other constants as given in the text, Section [V|. In all figures densities are in units 
where m = \^m Planck . 
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FIG. 2. logp = A and log(energydensity((^i)) = B versus the number of e-folds of expansion 
in reheating./? is the density of the radiation gas and tpi is the inhomogeneous (or particle) part of 
the inflaton field. The parameters are as for Figure 1. 
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FIG. 3. log(<| V 2 1)) = A and log(C) = B where tp is the metric perturbation and ( is defined in 
the text. The parameters are as for Figure 1. 
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FIG. 5. log(C) of Figure 3 shown in more detail for the first 8 efolds of reheating. 




3 4 
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FIG. 6. log(C) of Figure 4 shown in more detail for the first 7 efolds of reheating. Comparison 
of this with Figure 5 illustrates differences in ( in the cases with and without frictional decay of 
the scalar fields. 



APPENDIX A: BIANCHI IDENTITY AND THE FRICTIONAL DECAY TERMS. 



Working to first order in the metric perturbation ip a Bianchi identity gives 

D U T» = E P + V 'S V + X 'S X = 



(Al) 
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E p = ff + (p + p) (3a - 3ip' + {di - 2^ tl )u du i ) (A2) 

S v = (l-2i/;) (<p" + 3aip' - 3ip V) - tfV + a 2 m^ - (1 + 2^)p >M (A3) 

5 X = (1 - 2<4>)( X " + 3a X ' - 3^V) - + - (1 + 2^)x,i,< (A4) 

where, from the space-time Einstein equation, the fluid velocity, u l , is given by 

cT 2 (p + p)w M i = -<p V,i + (4ttG) _1 + aijj). ( A5) 

Equating the expressions in (|A2| , [A3| , [A~4"D to zero separately gives the fluid equation and the 



scalar equations. To add friction we use Eq. (|Al|) to write 

[E p - Y^'f - T x (x') 2 } + </[S v + I>'] + X '[S X + T xX '} = 0. (A6) 

This equation, which is equivalent to the Bianchi identity, is solved by equating each term 
in square brackets to zero. This modifies the fluid and scalar equations to give a transfer of 
fluid and scalar densities which however is modulated by the complexity of our complete set 
of equations. In our model p = po+Pi, <fi = ¥>o = <Pi where po, <po are homogeneous and pi, <pi 
are non-homogeneous and along with i/j linear in the same stochastic variables - this being 
forced through consistency of the various equations of the theory; x is linear in different 
stochastic variables. Each of the first two equations can be divided into a homogeneous 
part and a part linear in the stochastic variables with the occurence of quadratic or cubic 
stochastic variable terms being dealt with by appropriate ensemble averaging. Thus each of 
the first two equations is divided into two two halves which must separately be zero and after 
some cross substitution we get, together with the equation from the third square bracket 
the equations for po, pi, <fo, <fi, and x which appear in the text. 



APPENDIX B: EVALUATION OF THE LOOP INTEGRALS. 



Integrals such as that in Eq.(pO[) occur throughout the equations of motion and are 
evaluated numerically and we have to adopt a finite range of wave number, k. If the integrals 
diverge as k — > oo then the upper limit of the k integration forms a cut-off which is the 
crudest way of dealing with such ultra-violet divergences. Two points of view may be taken 
on this. Firstly this may be considered equivalent to a renormalization procedure in mass and 
other quantities. Secondly the Lagrangian used may be considered as an effective Lagrangian 
which has absorbed extra degrees of freedom coming from supersymmetry which eliminate 
divergences at higher momentum, the cut-off representing this effect. We have taken a cut-off 
to correspond to a wavelength of H^ 1 so that k cuto ff w 2naH. 

Also an infra-red cutoff is needed when the spectral index n is less than 1 irrespective 
of the particular inflation model (of the normal type); for the power law inflation model n 
is always less than 1. Inflationary theories which are in agreement with COBE observations 
have a value of n near, either greater or less than, 1 and thus there is a steep rise of 
(for a principal example) ipk as k — > as mentioned above in section [IV B 4|. First we 
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should make the general remark, independent of any particular value of n, that for k = 
the inflation era perturbation is homogeneous and may be regarded as incorporable in the 
classical homogeneous field <^o- The problem arises that for k near to the perturbation 
is nearly homogeneous and also a singularity of the loop integrand can arise at the end of 
inflation and in reheat. When the actual integral diverges the problem is acute and we deal 
with it by an infra-red cut-off taking the lower limit of the integration as k = ko so that Eq.[l8| 
becomes (ipf) = (2n)~ 3 J k d 3 k{pk<pl. The perturbations for k < k we regard as hidden by 
renormalization of the field (pQ. (In the equations of motion (cpf) always occurs in association 
with <Pq.) k should be small enough so that a (p^ which gives rise to observable effects is 
explicitly included in the loop integral. Thus k < k m = a m H m where the suffix m denotes 
the value at the time when the matter era begins, as observations of the CMBR fluctuations 
are possible at fc-values nearly of this order; k m is tiny compared with the corresponding 
value of aH at reheat. 

Now consider any value of k such that k <C aH throughout the reheat era; then 
(pk oc k n l 2 ~ 2 at the end of inflation and since for such small values of k (wavelength much 
greater than the Hubble radius) the k 2 term in the reheating equations of motion can be 
neglected means that the initial proportionality is maintained throughout reheat and thus 
cp k = (p{j)k n l 2 ~ 2 . We thus evaluate the loop integral as 

(tpl) = (27T 2 )- 1 f" 1 <p(T) V {T)*k n - 2 dk + (2tt)- 3 / a*k W l (Bl) 

J ko J k\ 

where k\ <C aH throughout reheat and the last term can be evaluated approximately using 
discrete values of the integrand since the delicate feature is included in the first term on the 
right hand side. 

We shall now treat these small k contributions in more detail, distinguishing the two 
cases n < 1 and n > 1. 

Case l:n < 1. Let n = 1 — e where e > 0. The analytic integration in Eq.jBl] is 

f kl k- l ~ e dk = (V - K e )/t. (B2) 

Jko 

Because k = 0{k m ) and e = O(.l) the first term on the right of Eq.[Bl] contributes a 
relatively large positive quantity to the loop integral. For the power-law inflation model 
n < 1 so this is the case that applies and is used in this paper. 

Case 2: n > 1. Let n = 1 + e where e > 0. 

( kl k- 1+e dk = {k\ - k e )/e w k\/e. (B3) 

Jko 

So the corresponding contribution is much smaller than that in the preceeding case which 
applies in this paper. 
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